Gravitational Lenses 
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Consider a photon moving past a point of 
mass, M, with an starting “impact param- 
eter,” €. From classical Newtonian grav- 
ity, the photon will undergo an acceleration £ 
perpendicular to the direction of its motion. 
Under the Born approximation, the amount 
of this deflection can be calculated simply 
by integration along the path. 
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If we substitute dx = cdt, then i 
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The integral in (21.02) is analytic, and works out to 2/€?. So 
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and the Newtonian deflection angle (in the limit of a small deflec- 
tion) is 
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In the general relativistic case, however, gravity affects both the 
spatial and time component of the photon’s path, so that the 
actual bending is twice this value. Thus, we define the angle of 
deflection, otherwise known as the Einstein angle, as 
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To understand the geometry of a gravitational lens, let’s first de- 
fine a few terms. Let 


distance from the observer to the lens 


Dy = 

Ds = distance from the observer to the light source 
Das = distance from the lens to the source 

8 = true angle between the lens and the source 

0 = observed angle between the lens and the source 
E = distance from the lens to a passing light ray 

Q = the Einstein angle of deflection 


Note that 6, 0, a, and € are all vectors, and calculations must 
deal with negative, as well as positive angles. For a simple point- 
source lens, the vectorial components of the angles make no differ- 
ence, but in systems where the lens is complex (i.e., several galax- 
ies/clusters along the line-of-sight) the deflection angles must be 
added vectorially. 


In practice, a, 6, and @ are all very small, so small angle approx- 
imations can be used. Also, Dg, Ds, and Das are all much bigger 
than €. Thus, the deflection can be considered instantaneous (i.e., 
we can use a “thin lens” approximation). 


The geometry of a simple grav- 
itational lens system is laid out 
in the figure. From AOSI and 
the law of sines, 
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Since all the angles are small, 
sin(0—68) ~ 0—8, and sin(180— 
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(Again, for single point lenses, 
the fact that the angles are vec- 
tors are irrelevant. For sim- 
plicity, I will therefore drop the 
vector signs for the rest of these 
derivations.) 





Since we have already seen that 
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and, from the simple geometry 
of small angles, 6 = €/Da, 
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In other words, we have a relation between 6 and 0. But note: 
for a given value of 8, there is more than one value of 0 that 
will satisfy the equation. This is a general theorum of lenses. 
For non-transparent lenses (such as the Schwarzschild lens being 
considered) there will always be an even number of images; for 
transparent lenses, there is always an odd number of images. 


To simplify, let’s define the characteristic bending angle, ao, as 
a quantity that depends only on the mass of the lens and the 
distances involved 
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If 8 and 0 are co-planar (as they will be for a point source lens), 
then (21.10) is equivalent to the quadratic equation 
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whose solution is 
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This defines the location of the gravitional lens images as a func- 
tion of the true position of the source with respect to the lens, 2, 
and ag. 


Magnification for a Schwarzschild Lens 
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The relation 
1 
0 = 5 (5+ y/103 + 2) 2112) 


implies that at least one image will be magnified. To see this, 
let’s lay out a polar coordinate system with the lens at the center, 
and consider the light passing through a differential area, dA. At 
the lens, this area element is dA = BAġ Ap. However, due to 
the gravitational lens, the angles are distorted, so that the area 
the observer will see this light squeezed into area dA’ = OA@d A9. 
Thus, the lens has “focussed” the light from area dA to area dA’. 
The magnification will therefore be the ratio of the two areas 
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The last term is just the derivative of (21.12), or 
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which means that the magnification is 
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Note the limits: as 6 —> oo, u+} — 1 and u- — 0; in other 
words, there is no lensing. On the other hand, when 8 — O, 
both u, and u— go to infinity. As you can see from the curves, 
when 8 < ao, gravitational lens magnification is important. 


Another feature of this magnification is that it only occurs be- 
cause the solid angle through which the light is flowing is being 
distorted. Moreover, nothing happens in the @ direction; only the 
0 coordinate is getting compressed. Consequently, the shape of 
the source will be distorted. The generalized version of (21.13) is 
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where wọ is the solid angle subtended by the source in the absense 
of lensing, and the amplification is given by the determinant of the 
Jacobian matrix. Note that this Jacobian can either be larger than 
one (magnification) or smaller than one (de-magnification). It can 
also be positive or negative; this gives the images orientation with 


respect to the original source. (This is also known as the parity.) 
A general theorum of transparent lenses is that there is always at 
least one image with positive parity and magnification. 


Finally, note that the determinant in the above equation can be 
zero. (Or, equivalently, the magnification given in (21.15) can 
be infinite. These locations are called caustics. In practice, the 
magnification never gets to be infinite because a) the background 
sources have finite size, and b) if it ever does get near infinite, 
interference effects will become important. 





There are a few other relations that should be kept in mind. Since 
the magnification of a lens is given by 
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the brightness ratio of two images is 
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which simplies through (21.12) to 
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Some other convenient relations (which can be proved with some 
simple algebra) are 
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and 
0102 = —a%, (21.21) 


(Again, remember that the angles are vector quantities, so angles 
such as 0 can be positive or negative.) 


General Scale Lengths for Lenses 


It’s useful to keep in mind some of the numbers that are typical 
of gravitational lens problems. Remember that the key value is 
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For the MACHO project, which tries to detect stars in the Large 
Magellanic Cloud (D, = 50 kpc) that are lensed by stellar-type 
objects in the Milky Way, 
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In this case, the image separations 6 ~ ag are sub-milliarcsec, so 
you will never be able to resolve the two separately lensed images. 
But, with both images landing on top of each other, the source 
can be (greatly) magnified. 





For quasar lensing by a typical large elliptical galaxy (102° Mo) 
in an Einstein-de Sitter universe, 
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where h is the Hubble Constant in units of 100 km/s/Mpc. Here, 
the typical scale-length is a little over an arcsec, so two separate 
images can be resolved. 


Gravitational Lens Time Delay 


If the source behind a gravitational lens is variable, then, of course, 
the images we see will also be variable. However, the images 
will not necessarily vary at the same time: in general, there will 
be a time delay between the two events. This comes from two 
components. 


First, there is a purely geometrical time delay. Consider the value 
of Das + Da— Ds for each image of the lens. From AOSI and the 
law of cosines 


D? = D + DŻ, — 2D4Das cos(180 — a) (21.24) 
If we substitute 
cos(180 — a) = —cosa = — {1 — 2sin* (a/2)} (21.25) 
then 
D? = DŻ + Di, + 2DaDas — 4DaDaz sin’ (a/2) (21.26) 
or 
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= (Da + Das — Ds) (Da + Das + Ds) (21.27) 
If we now use the law of sines, 


sina _ sin(@— p) (21.28) 
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and consider that all the angles are small, so that sina & a, 
sin(@ — B) ~ 0 — B, sin? a/2 ~ a?/4, and Da + Das + Ds & 2Ds, 
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then 
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Since there are two values of 6, there are two values for the geo- 
metrical time delay. So the delay of one image with respect to the 
other (which is all we can observe), is 
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If you now recall that 
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then, we can substitute these in to get 
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But, we also have the relation 
1 
d= 5 (a= yite) (21.12) 


0102 = -a4 (21.21) 


which gives 


So, with this substitution, we get the final expression for the ge- 
ometrical time delay from one image to the other. 
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The second source of signal time delay arises from the passage of 
the photons through a potential. Any time light passes through 
a gravitational potential, it experiences the Shapiro time delay 


At = = Joa (21.34) 





ae 
where dl is the photon’s path through the potential well. For the 


case of a weak potential and small deflection, the delay a photon 
experiences going from source S to point I is 
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Similarly, the delay in going from point J to the observer is 
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So the total time delay for the potential is 
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We can now calculate the difference in the time delay from one 
image to the other, by realizing that at very large distances from 
the lense, the potential delay is negigible. So, if we pick some 
distance, D’ >> £, but D’ < Das, Ds, then we can write (21.37) 
as 
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The second and fourth terms do not depend on € and are the same 


for both images. So, when we take the difference of the two time 
delays we get 
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If we then substitute for € using 0 = €/Dg and note that D’ is the 
same for both rays, we get the final expression for the difference 
between the two Shapiro delays 
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The total time delay of one image with respect to the other is 
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Note, however, that this equation neglects two important effects. 


The universe is expanding an possibly non-Euclidean; thus, the 
variables Ds, Ds, and Dg, are complex functions of z, Qo, and 
A. Specifically, Dg and D, are angular distances, and Dg; is the 
angular size distance of the source as seen by an observed at the 
redshift of the lens. 





Cosmological lenses are not point sources. The best lenses are el- 
liptical galaxies (which are fairly concentrated collections of stars), 
but these galaxies are located inside of clusters. So real lenses in- 
clude the effects of both the galaxy and the cluster (which, of 
course, add vectorially). 
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The Uses of Gravitational Lenses 





Gravitational lenses are extremely interesting for four reasons. 





Because the magnification properties of lenses depend on the mass 
of the lensing object, microlensing experiments, such as OGLE 
and MACHO (which monitor millions of stars in the Galactic 
Bulge and the Large Magellanic Cloud), can probe the properties 
of halo stars (and other objects) that are two faint to see. In 
fact, by studying the light curve of a gravitational microlens, one 
can deduce not only the mass of the intervening lens (with some 
assumption about the lens’ distance), but also whether the lens 
is a single star or something else. Microlensing experiments have 
measured light curves that must be due to binary-star lenses, and, 
in a few cases, have detected planets around stars. (In theory, 
the amplitude of the microlens event does not depend on mass, 
so even Earth-size (and lunar-size) bodies can be detected. In 
practice, while the amplitude of the event is independent of mass, 
the duration of the event is proportional to mass, so, in order to 
detect planet-sized bodies, excellent and continuous time coverage 
is needed. 








Because of their high magnification, gravitational lenses can be 
used to observe high-redshift galaxies that would normally be 
much too faint to observe otherwise. There is a slight problem 
with this, in that the image of the galaxy is usually distorted 
(spread out) over many pixels (into an arc). However, if the light 
can be collected, you can use it to study a representative sample 
of normal (not exceptionally bright) high-redshift galaxies. This 
property has been used to detect and measure galaxies as distant 
as z > 7. 
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Because the time delay of a gravitational lens depends on dis- 
tance, these objects can be used to measure the Hubble Constant. 
Unfortunately, this method has two major drawbacks. First, the 
time delay depends not only on the distance to the lens, but its 
mass as well. This mass is not known a priori: for a typical el- 
liptical galaxy lens, one must try to infer the lens mass from the 
velocity dispersion of its stars. Second, the geometry of the lens 
is usually complicated. Elliptical galaxies are in clusters, and the 
cluster’s mass will also contribute to the time delay. One must 
therefore measure the velocity dispersion of the galaxies in the 
cluster in order to determine the cluster’s mass (while, of course, 
hoping that the cluster is virialized, and that none of the objects 
you are measuring are foreground or background objects). More- 
over, for distant sources, it is likely that the calculations will be 
further complicated by the presence of multiple lenses along the 
line-of-sight. 














All clusters (and halos) along a line-of-sight cause some (small) 
amount of magnification. Moreover, because the mapping of the 
source into the image plane is asymmetric, each event will cause 
the background images to be slightly distorted. ‘Thus, it is possible 
to obtain a census of galaxies (and clusters) just be looking for 
small distortions in the shapes of distant galaxies. By seeing how 
this (weak) signal changes with position on the sky, one can derive 
information about galaxy clusters throughout the universe, and 
constrain the evolution of heirachical clustering (which depends 
sensitively on A). 
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General Comments on Real Lenses 
As we have seen, the bending angle for a lens is 
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This implies that as € —> 0, the bending angle (and magnifica- 
tion) becomes infinite. This doesn’t happen with cosmological 
lenses, because the lenses are never point masses. Instead, the 
mass is distributed, i.e., M(&), which goes to zero as € goes to 
zero. So the equations do not diverge. 





In practice, the gravitational deflection due to a galaxy’s cluster 
is almost as important as that of the galaxy itself. Consequently, 
one must model the lens with a galaxy plus cluster potential, and 
this breaks the symmetry of the problem. Hence, the equation to 


solve is Se ae ets 
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In other words, the bending angle is the vector sum of the gravi- 


tational deflections caused by all the mass in the two-dimensional 
cluster. 
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The best lenses are massive and compact objects. For cosmo- 
logical lenses, this means elliptical galaxies. (They are the most 
compact of all galaxies.) Consequently, even though the majority 
of galaxies are spirals, it is generally reasonable to assume that 
all gravitatational lenses are caused by ellipticals. Ideally, to use 
a gravitational lens as a cosmological probe, you would like 


Optical identifications for both the lens and the source. (Occa- 
sionally, the lens is too faint to be studied.) 
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A measured velocity dispersion for the lens (to estimate its mass). 
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A variable source (to measure the time delay for the two images). 
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An angular separation between the two images of 1” < A9 < 2”. 
If the separation is any smaller, it is hard to resolve the images; 
if the separation is any larger, then the cluster’s mass begins to 
dominate the problem. Also, recall that 

) ) (21.41) 
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hence the gravitational time delay AT x A6?. Small separations 

keep the time delay manageable. (For instance, the lens 0957+561 

has a separation of 6”, and a time delay of between 1 and 2 years.) 
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5) VLBI radio structure, so that the shape of the images can be used 
to constrain the mass distribution of the lens. 





6) Multiple images to constrain the geometry. Four images are better 
than two, and an Einstein ring is best. 
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Statistical Microlensing 


Consider a Schwartzschild (point source) lens at redshift z in front 
of a source at redshift z,. The magnification of the two images is 
given by 
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If the two images are unresolved (i.e., the case of microlensing), 
then the total magnification is 
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With a little bit of math, this equation can be inverted to give the 
angle 6, which yields a total magnification of ur 
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The cross-sectional area for a magnification of greater than ur to 
occur is therefore 
EUT; Z, Zs) = np? (21.45) 


This means that at redshift zs, any object located in an area of 
olur, Z, Zs) = TD? 8? (21.46) 


will be microlensed. 
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Now, let’s calculate the probability that a source at redshift z, will 
be microlensed. To begin, let’s assume that the density of lenses 
is given by 7. (Since the best lenses are elliptical galaxies, you 
consider 7 to be the density of elliptical galaxies.) For simplicity, 
let’s also assume that 7 does not evolve with time (there are as 
many elliptical galaxies today as there were at redshift z), and 
that each lens has identical lensing properties (same mass, same 
mass concentration, etc.) 








If no is the density of lenses today, then the density of lenses at 
redshift z is given by 


n(z) = no(1 +z)? (21.47) 





and the total number of lenses between z and z + dz is 
dN(z) = ndV = no(1+ 2z)?dV (21.48) 


where dV is the cosmological volume element. (For Euclidean 
space, this would normally just be 4rr?dr, but because of the 
universal expansion and the possible non-flat geometry of space, 
the actual expression is a complicated function of z, Qo, and A.) 


Now if the density of lenses is small enough so that the probability 
of two lenses overlapping each other is negligible, then the total 
area at redshift z, that is microlensed up to at least a magnifica- 
tion of ur is just the sum of the cross-sections of all the lenses, 
oe 
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The probability of a microlens is then just this area, normalized 
to the total surface area of the universe at redshift zs 


Ay = 40 Ding (21.50) 


where D2, , Is the angular-size distance of redshift z. 
This “statistical microlensing” is extremely important. In fact, 
with typical numbers for the space density of elliptical galaxy 
lenses, you find that by z = 2, half the objects in the sky are 
magnified by at least 10%! 
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